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It. Abatraat 

The overall objective of this effort In support of the Magsat 


project Is to study the feasibility of modeling magnetic fields due | 
to certain electrical currents flowing In the earth's Ionosphere and 
magnetosphere. This second quarterly status and technical progress 
report discusses efforts devoted to reading Magsat data tapes ''n 
preparation for further analysis of the Magsat data. The report 
also describes a new modeling procedure that Is being developed to 
compute the magnetic fields at satellite orbit due to hypothesized 
current distributions In the lonosp' ''re and magnetosphere. This | 
technique utl.^lzes a linear current element representation of the i 
large- G''.ale space-current system. Several examples of the model I 
field perturbations computed along hypothetical satellite r-rblts | 
are shown. I 
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I. IMTItODUCTION 


The overall goal of this investigation is to study the feasibility 
of nodellng magnetic fields produced by certain electrical currents 
flowing in the earth's ionosphere-magnetosphere systoi. Vector mag- 
netic field measurements from the near-polar oblting Magsat satellite 
contain, in addition to the main geomagnetic field and crustal anomaly 
fields, components that arise from these external currents. In meet- 
ing the ultimate goals of the Magsat project, it is desirable that the 
external current effects be identified in the observations, and sub- 
sequently separated from the internal field. The objective of this 
effort will be to determine the capability of a modeling procedure to 
facilitate the separation of these external and internal components. 

The approach of this feasibility study shall be to develop forward 
modeling procedures through which the magnetic effects of model currents 
may be derived. It is Intended to Include, separately, the equatorial 
electrojet, currents, and the effects due to auroral zone and polar 
cap currents including the high latitude ionosphere-magnetosphere 
coupling currents. In each case candidate current systems will be de- 
vised and resulting "typical" magnetic field signatures calculated for 
comparison with Magsat observations. The grid-cell modeling procedures 
developed by J.L. Klsabeth provide one method for calculating certain of 
these currents. This as well as ocher techniques will be evaluated. 

II. ACCOMPLISHMENTS DURING REPORTING PERIOD 

1. Data Reduction 

In Che first quarterly refort, we reported the development of soft- 
ware programs to read the Magsat Chronicle format data capes on a U.T.D. 


PDF 11/45 coaq>utcr, and printout aithar tha orbital data alona or both 
orbital and oagnatlc flald-valuaa fr<m both scalar and vactor awgnato- 
aaters for any spaclfled tlaa pariod contalnad on the sourca tape. 

The prograa cnaputes geodetic longitude and latitude, and altitude of 
the spacecraft and outputs this Information along with Inertial and 
magnetic coordinate positions. The magnetic "eld observations for 
each second are scanned with maximum, minimum, and average values for 
each scalar heat! and each vector component being printed. This soft- 
ware package Is now being extended so that selected portions of Magsat 
data can be stored on disc. Additional data reduction software will 
access the stored data, subtract a model field, and plot the resultant 
magnetic component deviations on a high resolution Interactive vector 
graphics terminal. This will permit direct compariseas between the 
measured magnetic field perturbations and the model results described below, 

2. Field Modeling 

A technique is being developed to model Che magnetic effects at 
Magsat altitudes due to hypothetical currents In and above the high 
latitude Ionosphere. In this forward modeling technique prototype 
current systems representing the auroral electrojet and Ionospheric 
polar cap currents as well as the field-aligned currents that link these 
low altitude currents to the distant magnetosphere are devised. This 
current system Is Chen represented for computational purposes as an 
array of linear current elements having finite length and a finite 
diameter with a pre-determlned cross-sectional current density distri- 
bution. For each hypothetical observation point Che three vector 
components of the resultant magnetic field are Co be computed as a 
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•uperposition of the contributions due to each current element In the 
vorld-space. The resultant magnetic perturbations for each vector- 
field component are displayed on a high resolution vector graphics 
terminal by means of a computer program designed to allow the operator 
to interactively modify the model paraiieters. The initial stages of 
development of this model dlscussec' in the first quarterly report had 
been restricted t** a .lypothetlcal satellite orbit at 90” Inclination 
in the dawn dusk meridian plane and the initial computations Included 
only the sunward component of the perturbation field at \'arlous alti- 
tiides. 

During this second quarter the capabilities of this model have 
been extended considerably. The model is now capable of computing 
three vector components of the magnetic field arising from an assumed 
current system. The restriction on satellite orbits has now been 
almost entirely eliminated allowing the magnetic field components to 
be computed for virtually any satellite orbit over a range of inclina- 
tions and altitudes and having an arbitrary angle between the orbital 
plane and the earth-sun line. Figures 1 through 4 show sample plots 
of the model magnetic field output for a satellite at 500 km 
altitude with various orbital inclination angles. The main plots on 
these figures show three orthogonal components of the magnetic field 
each cal'ulated at 94 separate observation points as a function of 
latitude. The solid curve is the sunward component, the dot-dash curve 
is the dawn-dust component and the dashed curve the vertical component. 
The input current system is a "classical" large-scale Blrkeland sheet 
current model with downward directed currents in the high latitude 
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postmldnlght and tha low latitude pre-midnight portions and upward 
directed currents In the high latitude pre-midnight and low latitude 
postmldnlght sectors. This current system is represented computation- 
ally In this example by 324 linear current elements as described 
above. The field-aligned currents are closed by M-S currents In the 
Ionosphere at 110 km altitude. No E-W ionospheric currents have been 
Included. The centers of the current sheets are located In the figures 
by the vertical lines at 61° and 70° Invariant latitude. The latitudinal 
distribution of current Intensity Is plotted along the top of the main 
panel. The clock dial In the lower right-hand quadrant depicts the 
satellite orbit. 

Two abstracts have recently been submitted for presentation at 
forthcoming scientific meetings. The purpose of these papers Is to 
describe the magnetic modeling technique developed under this contract 
and to present some of the Initial resulting magnetic perturbation 
signatures that have been derived from various model configurations 
of the polar Blrkeland and Ionospheric current systems. The first 
paper by Greer and Klumpar entitled "A Method of Calculating Magnetic 
Fields Due to Systems of Distributed Currents" has been submitted for 
presentation at the spring meeting of the American Geophyclcal Union 
and will primarily describe the modeling technique. The second paper 
entitled "A Technique for Modeling the Magnetic Perturbations produced 
by Field-Aligned Current Systems" by Klumpar and Greer will present 
results of several trial runs of the modeling program and attempt to 
show how these compare with the typical high latitude perturbations 
observed at Magsat. This latter paper will be delivered at the 
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Fourth Scientific Assembly of the Internstlonsl Associstlon of Geo- 
Bsgnetlsm and Aeronomy In August 1981. Copies of these two abstracts 
are appended to this quarterly report. 

III. PROBLEaffi ENCOUMTERED ANP REgamBTOATIONS 

No major problems have been encountered. 

IV. PLANS FOR NEXT REPORTING INTERVAL 

During the forthcoming quarter the development of data reduction 
software will continue. Programs to store selected portions of Hagsat 
data on disc file for subsequent analysis and plotting are now being 
written. In particular, routines will be developed and Integrated Into 
our existing software to compute the geopotential magnetic field at 
Magsat and to subtract away this model field. The resultant observed 
perturbations in the scalar magnitude and each of the 3 vector compon- 
ents will then be plotted by Interactive graphics routines designed to 
allow flexibility in plotting format. 

Further development of the linear element field modeling procedure 
will take place during the next reporting period. The emphasis will be 
placed upon Increasing the flexibility to choose more diverse Initial 
current configurations. 

A third major effort during the next quarter will be to obtain and 
bring into operation the grid cell modeling programs developed by J.L, 
Klsabeth. In the long term this grid cell modeling technique may be 
compared and contrasted to the linear elemental current model being 
developed at U.T.D. and discussed above. Ultimately some combinations 




of the two techniquM utilising the relative strengths of each might 
yield the isost realistic external field perturbetions in the Magsat data. 
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APPEHDICES 


Abstracts of papers suteittad for prasantatisn at scientific aeetings: 


A. A Method of Calculating Haaaetic Fields Due to Systems 

of Distributed Currents , suimltted tc Aaerican Geophysical 
Union, Spring Iteeting, BaltiMre, Md. Hay 2S>29, 1981 

B. A Technique for Hodelina the Masnetic Perturbations Pro - 
duced bT Field" Aliened Current Sv^ tews. sufaaltt^ for 
presentation at Fourth Scientific A^?«^ly of lACA, 
Edinburgh, Scotland, August 3-15, 1981 
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Appcndix A 


ABSTRACT 


SUBMITTAI. INFORMATIOH 


A METHOD OF CALCULATIHG MAGMETIC FIELDS 1. Spring Mtstltig 

DOE TO SYSTEMS OF DISTEIBDTED CnSBZHTS 


D. M. Gror 

D M. Elmpar (both at: Cantar for Space 

Sclancaa, Dolvarslty of Taxaa at Dallaa, 

Box 688, Rlchardaon, Taxaa 75080) 

Electrical currants in the ionoaphara and 
in the maffiatoaphare produce large amplitude 
magnetic field perturbations that are detect- 
ed by the highly sensitive magnetometers on 
the polar orbiting MAGSAT satellite. This 
paper describes a computationally fast and 
mathoatically simple method that has been 
developed and applied to modeling vhe 
magnetic field produced by the Blrkaland 
current system. The technique allows the 
magnetic field of a distributed current 
system to be calculated by representing such 
a system with an arbitrary number of hypo- 
thetical linear current elements. The 
facility of this method derives from the 
use of an analytical expression for the 
magnetic field of a straight current 
filament having an extended and smoothly 
varying current dedslty. The cross sectional 
current density profile of such a filament 
looks somcfdiat like a square wave pulse, of 
arbitrary width, with rounded comers. 

Thiis the system is free from unwanted 
discontinuities and the field component in 
any dlractlon and at any point in the model 
space Is easily calculated. Magnetic 
perturbations typical of thosa encountered 
at auroral latitudes by MAGSAT are produced 
by the model when realistic current configura- 
tions are chosen. 
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Appnulix B 


k TECHNIQUE FOR MODELING THE HA(^ETIC PERTURBATIONS PRODUCED 
BY FIELD-ALIGNED CURRENT SYSTEMS 

D.M. Klumpar «nd D.M. Greer (Center for Space Sciences, University 
of Texas at Dallas, Richardson, Texas 75080, U.S.A.) 

This paper presents results of a computational procedure that utilizes 
various assumed distributions of ionospheric and field-aligned currents to model 
magnetic perturbations observed at high latitudes frcmi the polar orbiting MAGSAT 
satellite. The highly sensitive vector magnetometers on MAGSAT repeatedly ob- 
serve magnetic field perturbations on essentially every transit of the high 
latitude ionosphere. These perturbations, with field conponents lying pre- 
dominantly in the magnetic East-West direction, are customarily vletMd as the 
signatures of oppositely directed paired sheets of electrical current flowing 
parallel to the geomagnetic field. These paired current sheets are typically 
regarded as being highly restricted in latitudinal extent and elongated in 
magnetic longitude. The model developed under this research utilizes a comput- 
ationally fast and mathematically simple technique that allows the magnetic field 
of a distributed current system to be calculated by representing stich a system by 
an arbitrary number of hypothetical linear current elements. The facility of 
the technique derives from the use of an analytic expression for the magnetic 
field of a linear current element having an extended and smoothly varying cross 
sectional current density; thus eliminating unwanted discontinuities. Mag- 
netic perturbations typical of those encountered at auroral latitudes by MAGSAT 
are produced by the model when realistic current configurations are chosen. 

Direct comparisons between the owdel field perturbations and those measured by 
the MAGSAT magnetometers permit more refined models of the Blrkeland currents to 
be developed. 
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